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Simplified model for current-voltage characteristics of weak links is sug- 
gested. It is based on approach which considers Andreev reflections as re- 
sponsible for the dissipative current through the metallic Josephson junction. 
The model allows to calculate current-voltage characteristics of weak links 
(superconductor - normal metal - superconductor junctions, microbridges, su- 
perconducting nanowires) for different thicknesses of the normal layer at dif- 
ferent temperatures. The current-voltage characteristics of tin microbridges 
at different temperatures were computed. 

1 Introduction 

Superconductor - normal metal - superconductor (SNS) junctions have the 
current- volt age characteristics (CVCs) with the rich peculiarities. Given cer- 
tain parameters of junction, CVCs of SNS junctions demonstrate the current 
peak at the small voltage, the excess current, the subgarmonic gap structure 
and the negative differential resistance at low bias voltage. Such nonlin- 
ear CVCs make SNS junctions to be promising for application to low- noise 
mixers in submillimetre-wave region PUI2]; switcher [3] or nanologic circuits 
111 

Description of CVCs of SNS junctions was subject of many articles and 
there were recognized key role of multiple Andreev reflections at NS interfaces 
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The main features of CVCs enumerated above are successfully described 
by Kummel - Gunsenheimer - Nicolsky theory (KGN) [7] where time depen- 
dent Bogoliubov - de Gennes equations are solved and wave packets of the 
nonequilibrium electrons and holes are considered. KGN theory is applica- 
ble for relatively thick and clean weak links where the normal metal layer 
N has the thickness 2a larger than the coherence length of superconductor 
and the inelastic mean free path I larger than 2a. Simplified model in frame 
of KGN theory was developed by L.A.A. Pereira and R. Nicolsky [TI]. This 
model is relevant for the weak links with thin superconducting banks S. The 
contribution of the scattering states [7j is omitted in [TT] . 

KGN and Pereira - Nicolsky model were applied earlier to describe the 
experimental CVCs [I2| E3 UH EH EH H7J- Experience of the applications 
demonstrates that oversimplified Pereira - Nicolsky model gives only qualita- 
tive description. New simple modification of KGN theory is proposed in this 
article. It is shown that CVCs of SNS junctions can be computed without the 
all complex Ansatz of KGN theory. I hope it will lead to more extensive using 
of the KGN based approach to the calculation of weak link characteristics 
than it was earlier (TJ E2 H d EES EE ETJ • 



2 Current -volt age characteristics 
2.1 Model 

Let us consider a voltage-biased SNS junction with a constant electric field 
which is in negative z direction perpendicular to the NS interfaces and exists 
in the N layer only (Fig. 1). The normal layer has the thickness 2a. The 
thickness of the superconducting bank is D — a » 2a. 

Accordingly KGN [7] the expression for CVC of SNS junction with thick 
superconducting banks (D — a » 2a) can be written as following: 
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Figure 1: SNS junction modelled. 

where m* is the effective mass of electron, g r (E) is the two dimensional 
density of states, Pn(E) is the probability of finding of the quasiparticles 
with the energy E in the N region of the thickness 2a, I is the inelastic mean 
free path and Rn is the resistance of the N region, A is the value of energy 
gap of superconductor at the temperature T, k z p is the z component of Fermi 
wave vector of quasiparticles, n is the number of Andreev reflections which 
quasiparticles undergo before they move out of the pair potential well. 

Eq.flTJ) is for the time averaged current that includes the voltage depen- 
dence only within the integral limits. 

2.2 Density of states 

To operate on Eq.(JT]) one should calculate the density of states [T5] : 

(2) 

where A is the normal layer area, k z p, r defines the value of k zF for which 
E r = E. 

The energy spectrum E r {k Z F) consists of the spatially quantized bound 
states and the quasicontinuum scattering states. The energy eigenvalue equa- 
tion for the spatially quantized bound Andreev states [7] is transcendental 
and calculated numerically only: 
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Er ^ = 2^F ( r7r + arccos ~£) ' (3) 

where r = 0,1,2,. . . 

Let us simplify Eq.©. The expansion of arccos(£' r /A) in © to Taylor 
series (tt/2 — E r /A +. . . ) up to second term and subsequent expressing of 
E r (k Z F) are executed: 




Figure 2: Energy of the bound Andreev state with r = 0; 2a = 5000 A; 
A = 0.57 meV; k F = 1.62 A" 1 . 1) Eq.©, C = 0; 2) Eq©, C = 1; 3) the 
exact solution of Eq.Q; 4) Eq.©, C = vr/2(l - am*A/h 2 k F ); 5) Eq.©, 
C = tt/2. 

Dependence E r {k zF ) © (curve 2) and the numerical solution of Eq.© 
(curve 3) are shown in Fig. 2. The better agreement with the numerical 
solution of Eq.© is attained by insertion of the correcting multiplier C 
before h 2 k zF /2am*A in Eq.©: 

Br(M ^ M)/ ( 1+c |^) (8) 
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If C = then the spectrum of Pereira - Nicolsky model is reproduced 
(curve 1, Fig. 2). R. Kummel used Eq.(j3J) with C = tt/2 [TU] for approx- 
imated calculation of the energy spectrum (curve 5, Fig. 2). I suggest the 
variable multiplier C = n / 2(1 — am* A/h 2 kp) for C > 1 and C = 1 otherwise. 
Such choice of C provides a good agreement of Eq.fjSJ) (curve 4, Fig. 2) with 
the numerical solution of Eq.(j3J) for different relation of a, m*, A, kp. 

The density of the bound states with © becomes 




For quasiparticles from the quasicontinuum states the energy spectrum 
is approximated by the continuous BCS spectrum of a homogeneous super- 
conductor [7J dH] : 
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Figure 3: Density of states g{E) of an SNS junction with thick supercon- 
ducting banks resulted by [18 j (dotted line) and g(E) calculated by Eq.© 
with C = 7r/2(l - am*A/h 2 k F ) (solid line). D = 70000 A; 2a = 5000 A; 
T c = 3.77 K; k F = 1.62 A^ 1 . 



5 



In the case of SNS junction with thick superconducting banks the effective 
energy gap A* equals A. The density of the quasicontinuum scattering states 
is 

g(E) = — ^ — Ty-kpD - 8 
yv ; 7T 2 h 2 VE 2 - A 2 V ; 

The density of states resulted is shown in Fig. 3. 
2.3 Current density 

The probability Pn(E) of finding of the quasiparticles with the energy E in 
the N region is given by Eq.(2.19) of [Zj: 

with the penetration depth A = -/jfjrp for E < A, \ < D — a and 
A = D — a otherwise. 

For the quasiparticles from the scattering states Pn(E) = 2a/ 2D. Let us 
accept for the sake of simplicity A >> a and therefore Pn{E) = 2a/ '2 A for 
the bound states. 

The current density of quasiparticles from the bound states is resulted 
with ©: 

. < em* 2 a 2 ^ f 2a 
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With using (jSJ) the current density of quasiparticles from the quasicon- 
tinuum states is resulted: 



dE V F r± == tanh(—=— ) (11) 
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Here and further E\ = — A + neV for — A + neV > A and E\ = A 
otherwise. 

Neglecting the small term I have 



e / 2a \ A+ f V Et&nh(E/2k B T) 
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If eV >> fcsT, A the integral in (|T2"j) can be transformed and the excess 
current density is resulted: 

fo (V^^A t a„h(^)exp(-f) (13, 

This excess current density is the same as one in KGN (Eq.(4.12) in [7j). 

Note that jbsiY) dependence does not change practically if the second 
summation in (jlOJ) is interrupted at r = 0. Therefore the total current 
density is 
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C = tt/2(1 - am*A/h 2 k F ) for C > 1 and C = 1 otherwise; E\ = —A + neV 
for —A + nel^ > A and E\ — A otherwise. 

Eq.([14|) is the main result of this simplified model. The model allows to 
calculate CVCs of weak links for different thicknesses of the normal layer at 
different temperatures. The subgarmonic gap structure, the excess current 
and the current peak at the small voltage are reproduced on CVCs. 

2.4 Comparison with experimental current-voltage char- 
acteristics 

The model was used to compute two sets of CVCs of tin microbridges. These 
detailed measurements of the current biased CVCs were performed by V.N. 
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Gubankov, V.P. Kosheletz, G.A. Ovsyannikov in 1977-1981 [2011211 (Fig. 4) 
and M. Octavio, W.J. Skocpol, M. Tinkham in 1978 22^ (Fig. 5). Both sets 
of CVCs have the similar peculiarities: the subgarmonic gap structure, the 
current peak at the small voltage and the excess current. 




Figure 4: The current- voltage characteristic of Sn microbridges. Experiment 
[2*T] (points) and calculations (solid lines). 

Comparison of the computed /(V) curves and the experimental V vs. / 
dependencies displays satisfactory agreement at temperatures smaller than 
0.99T C . Presented in Fig. 4 and Fig. 5 curves were calculated with the rea- 
sonable parameters: the critical temperature, the energy gap at zero temper- 
ature, the Fermi wave vector of Sn (T c = 3.77 K, A = 0.57 meV, kp = 1.62 
A -1 ). The length of microbridges 2a is 5000 A and I = 15a. The BCS 
dependence of A on T was used. 

The high voltages regions on experimental CVCs are close by the com- 
puted curves at higher temperatures. It is possibly reasoned by selfheating 
occurred at high voltages in these experiments |2*Tll22j . Some discrepancy of 
the computed curves and the experimental points at low voltages is because 
there were the current-biased CVCs in experiments instead voltage-biased 
one. Agreement of the model and the experimental CVCs disappears at 
temperatures near T c : the calculated current peak and the excess current is 
smaller than corresponding experimental currents (e.g. CVCs at 3.751 K in 
Fig. 4 and at 3.746 K in Fig.5). 
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Figure 5: The current-voltage characteristic of Sn microbridges. Experiment 
[22] (points) and calculations (solid lines). 

3 Conclusion 

Simplified model for calculation of current- voltage characteristics of the weak 
links (SNS junctions, microbridges, superconducting nanowires) was devel- 
oped. This model makes the KGN approach [7] more convenient for de- 
scription of experiments. The model was applied for computation of the 
current- volt age characteristics of tin microbridges at different temperatures. 
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